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Abstract. We present new determinations of the blue to red supergiant ratio (B/R) in young open clusters 
at various metallicities. For this purpose, we examine the HR diagrams of 45 clusters in the Galaxy and of 
4 clusters in the Magellanic Clouds. The identification of supergiants is based on spectroscopic measurements 
(with photometric counts to check the results). The new counts confirm the increase of the B/R ratio when the 

metallicity increases with the following normalized relation : nf/g; = 0.05 -e z ® , where Zq — 0.02 and (B/R)q 
is the value of B/R at Zq which depends on the definition of B and R and on the age interval considered (e.g. 
for spectroscopic counts including clusters with log ape between 6.8 and 7.5, (B/R)q = 3 when B includes O, B 
and A supergiants). 
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1. Introduction 

The variation with metallicity Z of the number of blue 
and red supergiants is important in relation to the na- 
ture of the supernova progenitors in different environ- 
ments (Langer 1991ab) and its effects on the luminous 
star populations in galax ies (e.g. Cervino & Mas-Hesse 
1999] ). This ratio also constitutes an 



1994; Origlia et al. 



important and sensitive test for stellar evolution models, 
because it is very sensitive to mass loss, convection and 
mixing processes (Langer & Maeder 1995 ). Stellar evo- 
lution models can usually, by adjustement of parameters 
such as the rate of mass loss by stellar winds, reproduce 
the observed ratio for a given metallicity, but are not able 
to reproduce its variation with the metallicity. Thus, the 
problem of the blue to red supergiant ratio (B/R ratio) 
remains one of the most severe problems in stellar evolu- 
tion. 

Since the last studies of the B/R ratio in the Galaxy 
and in the Magellanic Clouds more than fifteen years ago 
(Meylan & Maeder |1982| ; Humphreys & McElroy |1984| ), 
many new spectroscopic and photometric measurements 
have been performed. Our aim here is to reexamine this 
question accounting for these recent observational im- 
provements. We chose to use stellar clusters instead of field 
stars for obvious reasons: stars are at the same distance, 
with the same age and above all have the same chemical 



composition. Moreover, the knowledge of the age of the 
clusters enables us to estimate the initial masses of the 
supergiants used to derive the B/R ratio. 

Thus, the present work is in the continuity of the work 
of Meylan & Maeder ( |1982| ). However, the present study 
presents two significant improvements. First, the identifi- 
cation of supergiants is based on spectroscopic measure- 
ments instead of on photometric colours. Spectroscopic 
measurements allow a more accurate differentiation be- 
tween supergiants and main sequence stars (see however 
the discussion in Sect. 4). Second, the number of stellar 
clusters in our galactic sample is increased by a factor of 
four. 

In Sect. 2, we review the various studies about the 
B/R ratio. Section 3 presents the new counts. Our results 
are discussed in Sect. 4, and conclusions are given in Sect. 
5. 



2. Previous studies 



Walker (|196J) observed that the B/R ratio in M33 de- 
creases when the distance to the centre of the galaxy in- 
creases. Van den Bergh (1968) explained this variation 
by a radial metallicity gradient in the disk of M33, the 
metallicity being higher in the inner parts. Humphreys 
& Sandage (1980) reexamined the B/R ratio in M33 
and confirmed its decrease when the galactocentric dis- 
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tance increases. Freedman (1985) contested the results of 
Humphreys & Sandage (198C) concerning the variation 
of the B/R ratio in M33 and concluded that this was 
just an effect of incompleteness. However, Ivanov (1998) 
concluded that the diminution was real and was in good 
agreement with the radial metallicity gradient observed in 
M33 (McCarthy et al. |l995| ). 

Hartwick (|l970| ) noticed the decrease of the B/R ratio 
in the Milk y Way as th e galactocentric radius incr eases . 
Robertson (|l973| , |l974| ), Hagen & van den Bergh ( |l974|) 
compared theoretical evolutionary tracks to observed HR 
diagrams (HRD) of stellar clusters in the Galaxy and in 
the Large Magellanic Cloud (LMC). They attributed the 
observed differences between the HRD to differences in 
the metal content. 



3. Results from new cluster studies in the Galaxy, 
LMC and SMC 

3.1. Spectroscopic counts 



Humphreys (1979a) investigated the variation of the 
B/R ratio in the LMC. She concluded that the radial vari- 
ation of this ratio must be weak or absent, in agreement 
with the study of the chemical composition of the HII 
regions made by Pagel et al. (1978). In a study of the su- 



pergiants in the Galaxy and in the LMC, Humphreys & 



Davidson (1979) found a decrease of the B/R ratio when 
the metallicity decreases. Cowley et al. (1979) examined 
the variation of the B/R ratio across the face of the LMC 
and found that it increases by a factor of 1.8 when the 
metallicity increases by a factor of 1.2. 

The first study of the B/R ratio using stellar clusters 
in the Galaxy, the LMC and the SMC ( Small Magellanic 
Cloud) was made by Meylan & Maeder ( 1982| ), who con- 
cluded to a steep increase of B /R when the metallicity in- 
creases, with a difference of about an order of magnitude 
in B/R between the Galaxy and the SMC. Humphreys & 
McElroy (1984) reexamined the B/R ratio in the Galaxy, 
the LMC and the SMC, and gave new values which ac- 
counted for the incompleteness of the data. They found a 
B/R ratio about ten times higher in the central parts of 
the Galaxy than in the SMC. 

Langer & Maeder (1995) summarized some results 
about the values of the B/R ratio and gave a B/R ra- 
tio approximately six times larger in the solar neighbor- 
hood than in the SMC. After comparison of different stel- 
lar models with observations, they concluded that most 
massive star models have problems reproducing the de- 
crease of B/R when the metallicity decreases. Deng et al. 
(1996) used the su pergiant stars catalogue from Blaha & 
Humphreys ( 1989 ) to make new counts in the Galaxy and 
in the LMC. They found a B/R ratio about five times 
higher in the Galaxy than in the LMC. 

To summarize the B/R ratio appears to be an increas- 
ing function of metallicity. This trend is observed by stud- 
ies in the Galaxy, the Magellanic Clouds as well as in the 
Triangulum galaxy. In this work we shall take the oppor- 
tunity of recent improvements brought to our knowledge 
of Galactic and Magellanic clusters to reexamine carefully 
this question. 



We select in the Webda data base (Mermilliod 1995 ; 
irttp : / /obs www . unige . ch / webda ) young clusters with 
log age between 6.8 and 7.5, corresponding to masses at 
the turn off between 8 Mq and 24 Mq and to initial masses 
of th e sup ergiants between about 8 and 30 Mq (Meynet 
et al. 1993 ). Let us note that this interval of masses for the 
supergiants probes a region of the HR diagram below the 
observed upper limit for the luminosity of the red super- 
giants, i.e. below Mboi ~ —9.5 (Humphreys & Davidson 



1979, 1984) 



We only keep the clusters with spectroscopic measure- 
ments for the brightest stars in order to distinguish su- 
pergiants from main sequence stars. In this way, 45 clus- 
ters are selected. For each cluster, we count the number 
of supergiant stars of spectral type O, B, A, K and M. 
In a few cases, the spectroscopic class for the same star 
is different depending on the authors. In these cases, we 
adopt the most recent determination. The galactocentric 
distance R of each cluster is calculated from the values 
of its galactic coordinates and from its distance to the 
sun given by the Webda database. For each galactocen- 
tric radius, a metallicity is associated by taking a value of 
Z = 0.02 at the solar position (we take R Q = 8.5 kpc for 
the solar galactocentric distance) and an averaged radial 
metallicity gradient of —0.08 dex kpc -1 . This value has 
been chosen from various studies of the galactic gradient 
determined by means of spectroscopic and photometric 
indices of stellar populations in open clusters (see the dis- 
cussion by Alibes et al. 2001). 

The results are given in Table [l], where B represents 
the number of O, B and A supergiants, while R is the 
number of K and M supergiants. No, Nb, Na, Nk and 
Nm indicate the number of supergiants with spectral type 
O, B, A, K and M respectively. Each cluster age is given 
in logarithm and taken from the Webda database. The 
region noted GC refers to the region inside the solar cir- 
cle and includes all the clusters with a galactocentric dis- 
tance shorter than 8.5 kpc. Likewise, the GAC region cor- 
responds to the region outside the solar circle and includes 
the clusters with a galactocentric distance higher than 8.5 
kpc. 

On the basis of these results, we examined the vari- 
ation of the B/R ratio with the galactocentric radius by 
grouping open clusters in different distance bins. In Fig. 
two different binnings are shown. In the first one, clusters 
are grouped in three main galactocentric distance inter- 
vals (6.5-8.0, 8.0-9.0 and 9.0-11.5 kpc), in order to have 
the same number of clusters (15) in each bin. We also sep- 
arate the clusters in five distance intervals (6.5-7.5, 7.5- 
8.0, 8.0-8.5, 8.5-10.0 and 10.0-11.5 kpc) to check whether 
the variation of the B / R ratio observed when considering 
three intervals is still present (due to the small number of 
supergiants counted for each clusters, it is not pertinent 
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Table 1. Number of supergiant stars in galactic open clusters with log age between 6.8 and 7.5. The quantity R is the 
galactocentric radius and Z is the metallicity. GC means clusters toward the galactic centre, GAC toward the galactic 
anticentre. The identification of supergiants is based on spectroscopic measurements. B includes O, B, A supergiants 
and R includes K, M supergiants (see text). 
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to adopt a finer binning). The metallicity of each bin is 
obtained by averaging the metallicities of all the clusters 
in this bin. We precise here that the metallicities are not 
obtained from stars in each clusters, but from an adopted 
metallicity gradient of —0.08 dex kpc -1 . In this way, we 
assign a metallicity for each distance interval mentioned 
above (e.g. for the GC and GAC region we take respec- 
tively Z = 0.023 and Z = 0.016). The increase of the B/R 
ratio with the metallicity is shown in Fig. [j] and in Table ||. 

In order to obtain values of the B/R ratio for lower 
metallicities, we select young clusters in the Mag ellani c 
Clouds. According to a discussion by Maeder et al. ( 1999|) , 
the average value of the metallicity is Z = 0.007 for the 
LMC and Z = 0.002 for the SMC. Without sufficient 
spectroscopic data for the young clusters of the LMC, it 



is not possible to give a reliable value for B/R in the 
LMC. In the SMC, the only cluster satisfying the age cri- 
teria and having spectroscopic measurements is NGC 330 



(log age = 7.00 according to Cassatella et al. 1996). This 
cluster has been well studied and so different values ex- 
ist for the B/R ratio : B/R 
B/R = 9/15 (Carney et al. 



= 7/1 2 (Cayrel et al. |1988|) , 
198£ : Brocato & Castellani 



|1992| ; B omans & Grebel |1994D , B/R = 12/15 (Grebel & 
Richtler 1992 ). In these ratios, B only includes B super- 
giants. Thus, the value of the B/R ratio for NGC 330 lies 
somewhere between 0.5 and 0.8 (when only B supergiants 
are counted). The value of 0.6 is currently accepted and 
is reported in Fig. |l|. 
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Table 2. Number counts down to a limit magnitude M v u m . The values of (m v — M v )q and E(B — V) mentioned for 
the galactic clusters are from the Webda database (Mermilliod 1995| ) and from van den Bergh (1998) for the clusters 
in the Magellanic Clouds. B includes stars with (B — V)q < and R stars with (B — V)o > 1.2 (see text). 
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Table 3. B/R ratio for the galactic centre (GC), the galactic anticentre (GAC), SMC and LMC. B in the column 
called 'Spectroscopy' only includes B supergiants, while R includes K and M supergiants. In the other columns, B 
includes stars with (B — V)q < and R stars with (B — V)q > 1.2 (see text). 







Spectroscopy 






-3.25 


M vl 




-3.0 






-2.5 


M vl 




-2.0 




Z 


Bt 


R 


Bi/R 


B 


R 


B/R 


B 


R 


B/R 


B 


R 


B/R 


B 


R 


B/R 


GC 


0.023 


41 


10 


4.1 


189 


7 


27.0 


221 


8 


27.6 


308 


9 


34.2 


416 


13 


32 


GAC 


0.016 


29 


23 


1.3 


168 


17 


9.9 


206 


17 


12.1 


324 


21 


15.4 


460 


23 


20.0 


LMC 


0.007 








60 


24 


2.5 


79 


25 


3.2 


123 


27 


4.6 


165 


29 


5.7 


SMC 


0.002 


9 


15 


0.6 


67 


20 


3.4 


79 


21 


3.8 


139 


23 


6.0 


211 


40 


5.3 



f counting only B supergiants 



3.2. Photometric counts 

In order to complete counts based on spectroscopic mea- 
surements, we perform counts based on colour indices in 
the same way as Meylan & Maeder (1982). To this pur- 
pose, we select among the 45 open clusters mentioned in 
Table [l], those having photometric measurements down 
to M v = —2 and with at least 15 stars brighter than this 
magnitude. Hence 23 clusters are selected. The distinction 
between blue and red supergiants is made in the following 
way : when (B — V)o < 0, the star is considered to have 
a O, B or A spectral type and is included in B, whereas 



for (B — V)o > 1.2, the star is considered to be a K or 
a M and is included in R. This distinction remains valid 
for the SMC, where late- type supergiants are of earlier 
type than in the Galaxy and LMC (Humphreys 1979b; 
Elias et al. 1985). We then count all the stars in these two 
colour intervals down to a certain limit M v u m of the ab- 
solute magnitude. In order to be sure that the results do 
not depend on the value of this arbitrary chosen limit, the 
counts are made for four different limits : M v u m — —3.25, 
M vUm = -3.0, M vUm = -2.5 and M vUm = -2.0. The 
23 selected clusters and the results are listed in Table ^[ 
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Fig. 1. B/R ratio in the Galaxy and the SMC for clusters 
with log aye between 6.8 and 7.5. The distinction between 
blue and red supergiants is based on spectroscopic mea- 
surements. The triangles refer to B including O, B and A 
supergiants and correspond to five different distance inter- 
vals in the Galaxy (6.5-7.5, 7.5-8.0, 8.0-8.5, 8.5-10.0 and 
10.0-11.5 kpc). The dots refer to B including only B su- 
pergiants; they correspond to one value for the SMC and 
three distance intervals in the Galaxy (6.5-8.0, 8.0-9.0 and 
9.0-11.5 kpc). The dashed curve corresponds to the fit for 
B including O, B and A supergiants with (B/R)q = 3.0 
(see text). 



Because of the small number of clusters selected in this 
case, it is not pertinent to subdivide the clusters in sev- 
eral distance intervals or age intervals, and so we consider 
all the clusters with log age between 6.8 and 7.5 and only 
the two main regions GC and GAC to calculate the B/R 
ratio (see Table ||). 

For the Magellanic Clouds, we use the data from Elson 



(|1991|) , except for the SMC cluster NGC 330 (Va llenar i 
1994 ). Cluster ages are taken from Cassatella et al. ( |1990| ). 



The same criteria of selection as those used for galactic 
clusters give a total of 3 clusters selected in the LMC. 
In the SMC, NGC 330 is the only cluster satisfying both 
the criteria for age and for the number of stars brighter 
than M v = —2. The counts are made in exactly the same 
way as for the Milky Way. The selected clusters and the 
results of the counts are given in Table §, while Table | 
summarizes the values of the B / R ratio found in the four 
main regions (GC, GAC, LMC and SMC). Let us note here 
that our selected clusters have a mean age approximately 
the same in the GC, GAC and SMC regions (respectively 
7.0, 7.1 and 7.0) which is a nice feature for the comparisons 
performed here. 



Fig. 2. B/R ratio in the Galaxy for different age intervals, 
with distinction between blue and red supergiants based 
on spectroscopic measurements. B includes O, B and A 
supergiants. As in Fig. |l|, triangles correspond to five dif- 
ferent distance intervals in the Galaxy (6.5-7.5, 7.5-8.0, 
8.0-8.5, 8.5-10.0 and 10.0-11.5 kpc). Filled squares and 
circles correspond to three distance intervals in the Galaxy 
(6.5-8.0, 8.0-9.0 and 9.0-11.5 kpc) and open squares and 
circles correspond to the two main distance intervals GAC 
and GC. Triangles include clusters with log age between 
6.8 and 7.5 (i.e initial masses of the supergiants between ~ 
8 and 30 Mq). Filled circles refer to clusters with log age 
between 6.8 and 7.2 (initial masses of the supergiants be- 
tween 12-30 Mq), while filled squares include clusters with 
log age between 7.0 and 7.4 (initial masses of the super- 
giants between 9-18 Mq). Open squares correspond to 
clusters with log age between 6.9 and 7.1 (initial masses 
of the supergiants between 14-23 Mq) and open circles to 
clusters with log age between 7.1 and 7.3 (initial masses 
of the supergiants between 10-15 Mq). The dashed curve 
corresponds to the same fit as in Fig. |l|. 



From the results for spectroscopic and photometric 
counts, we find the following average relation : 



B / R ~nn C 3^- 

= 0.05 • e z o 



{B/R), 

where Zq = 0.02 and (B/R)q is the value of B/R at 
Zq which depends on the definition of B and R and on 
the age interval considered. For clusters with log age be- 
tween 6.8 and 7.5 (initial masses of the supergiants be- 
tween ~ 8 and 30 Mq), (B/R)q = 3.0 for spectroscopic 
counts when B includes O, B and A supergiants, and 
(B/R)q = 2.5 when only B supergiants are counted. For 
photometric counts, (B/R)q equals 17.2, 18.5, 23.0 and 
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24.0 for M vUm = -3.25, M vUm = -3.0, M vUm = -2.5 
and M v ii m = —2.0, respectively. 

4. Discussion of the results 

Fig. [I] clearly shows that the present study, based on 
the most recent spectroscopic measurements found in the 
literature, confirms the trend already found in previous 
works, namely the decrease of the B/R ratio with the 
metallicity. In the Galaxy, B/R is approximately 3 times 
higher in the GC than in the GAC (Table |). This vari- 
ation can also be observed by separating the Galaxy in 
three or five galactocentric distance intervals. Moreover, 
the B/R ratio in the SMC confirms the results obtained 
for the Galaxy, with a value approximately 2 times lower 
than in the GAC region. Fig. [l] shows also that the de- 
crease is stronger at high metallicity than at low metal- 
licity. T his is in good agreement with recent results from 
Massey ( 2002) ) who found a B/R ratio only slightly higher 
in the LMC than in the SMC. 

Concerning photometric counts, we can see in Table | 
that B/R is higher in the GC region than in the GAC 
region. The increase of the B/R ratio when the limit mag- 
nitude increases is simply due to the fact that more main 
sequence stars are included in B. The results are stable 
for the limit magnitudes —3.25, —3.0, and —2.5, for which 
B/R is approximately 2.5 times higher in the GC than 
in the GAC (in good agreement with Meylan & Maeder 



1982 and with the results from spectroscopic counts). For 



the limit magnitude of -2.0, the value of B/R for the GC 
seems relatively low, which likely reflects problems of com- 
pleteness. 

The B/R ratios in the LMC (Table |) confirm the 
decrease of B/R when the metallicity decreases. However, 
the results for the SMC seem less clear. This is probably 
due to the fact that ratios for the SMC are based on one 
single cluster (NGC 330). This is obviously insufficient to 
obtain precise B/R ratios and thus, in this case, results 
must be regarded with circumspection. 

What are the effects of grouping clusters of different 
ages on the results ? Let us recall that we selected clus- 
ters with log age between 6.8 and 7.5. This means that 
the B/R ratio we obtain are for supergiants with initial 
masses between about 8 and 30 Mq. Adopting narrower 
age intervals would give indication on the B/R ratio in a 
smaller range of initial masses. One can wonder to which 
extent the decrease of the B/R ratio with the metallic- 
ity would remain the same, if smaller age (mass) intervals 
were adopted. To check this point, we considered differ- 
ent age intervals (i.e. different intervals of initial masses 
for the supergiants). It is obvious that a finer age interval 
means less clusters selected and so a less reliable statistic. 
Therefore, we choose larger distance intervals when finer 
age intervals are considered. First, we calculate the B/R 
ratio in the same three distance intervals considered above 
(6.5-8.0, 8.0-9.0 and 9.0-11.5 kpc) for two age intervals : 
log age between 6.8 and 7.2 (initial masses of the super- 
giants between 12-30 -Mq) and log age between 7.0 and 



7.4 (initial masses of the supergiants between 9-18 Mq). 
We also choose two finer age intervals, log age between 
6.9 and 7.1 (initial masses of the supergiants between 14- 
23 Mq) and log age between 7.1 and 7.3 (initial masses 
of the supergiants between 10-15 Mq), for which we cal- 
culate the B/R ratio in the two main distance intervals 
GC and GAC. The results of these counts are shown in 
Fig. ^. We precise that the value of the galactocentric ra- 
dius R of each distance interval is obtained by averaging 
the galacocentric radii of the clusters found in the bin. 

We observe that for different age intervals, the B/R 
ratio always remains higher for larger metallicities. Fig. 
shows clearly that the variation of the B/R ratio with the 
galactocentric radius remains more or less the same what- 
ever the interval of ages (masses) considered. Moreover, we 
can see that considering finer age intervals changes only 
slightly the values of the B/R ratio. In that respect, the 
age interval between 6.8 and 7.5 appears to be a good 
compromise between the necessity to select a sufficiently 
high numbers of clusters to have a reliable statistics, and 
the necessity to somewhat restrain the domain of masses. 

Finally, let us stress that when comparisons with stel- 
lar evolution models are made, the same definitions of the 
blue and the red supergiants used for obtaining the ob- 
served B/R ratios has to be used to determine the theoret- 
ical B/R ratios. In that respect, it is worthwhile to recall 
here that even if a star is classified as a blue supergiant, 
it may belong to the upper end of the Main Sequence. 

5. Conclusions 

The decrease of the B/R ratio with the metallicity does 
not appear to be an artifact due for instance to incom- 
pleteness but instead a robust feature that stellar evolu- 
tion models should be able to reproduce. The fact that the 
present day grids of stellar models are still unable to ac- 
count for this feature demands some caution when evolu- 
tionary population synthesis models are used to interpret 
the integrated luminosity of high redshift gala xies. 

Recent studies by Maeder & Meynet ( |2001 ) show that 
the inclusion in the stellar models of the effects of rotation 
changes the B/R ratios. Physically, this results from the 
mild mixing, which leads to more helium in the region of 
the H-shell burning. The opacity is lower and the inter- 
mediate convective zone less important or absent. As con- 
vection implies a polytropic index n — 1.5, which means 
a relative compactness of the internal convective regions, 
the absence of an intermediate convective zone is a nec- 
essary condition for stellar expansion to the supergiant 
stage. Interestingly, rotating models are able to account 
for the numerous red supergiants seen at low metallicity, 
a feature that standard models could not reproduce. 
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